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Abstract 

An exact discontinuity capturing central solver developed recently, named MOVERS (Method of Optimal Viscosity 
for Enhanced Resolution of Shocks) m, is analyzed and improved further to make it entropy stable. MOVERS, 
which is designed to capture steady shocks and contact discontinuities exactly by enforcing the Rankine-Hugoniot 
jump condition directly in the discretization process, is a low diffusive algorithm in a simple central discretization 
framework, free of complicated Riemann solvers and flux splittings. However, this algorithm needs an entropy fix to 
avoid nonsmoothness in the expansion regions. The entropy conservation equation is used as a guideline to introduce 
an optimal numerical diffusion in the smooth regions and a limiter based switchover is introduced for numerical 
diffusion based on jump conditions at the large gradients. The resulting new scheme is entropy stable, accurate and 
captures steady discontinuities exactly while avoiding an entropy fix. 
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1. Introduction 

Of all the numerical methods developed to solve the hyperbolic systems of conservation laws of fluid dynamics 
numerically, approximate Riemann solvers have been the most popular. However, these approximate Riemann solvers 
are discovered to introduce shock instabilities like carbuncle shocks, odd-even decoupling, kinked Mach stems in the 
solutions 12] and they typically require addition of numerical diffusion in the form of entropy fixes to make them 
entropy stable. It is often nontrivial to extend these approximate Riemann solvers to hyperbolic systems other than the 
Euler equations of gas dynamics and further they are heavily dependent on the eigen-structure of the hyperbolic sys¬ 
tems. This dependency is due to the enforcement of strict upwinding procedure which usually results in sophisticated 
procedures of Riemann solvers or complicated flux splittings and are dependent on the projection of the solution on 
to the space of eigenvectors and characteistic decomposition. In this context, in recent research, attention is focused 
on developing simpler central solvers which are competitive to upwind schemes in terms of accuracy and efficiency 
of capturing essential flow features like shock waves, contact discontinuities and expansion waves E a 13 [Dig]. 

Of the new central schemes, the Method of Optimal Viscosity for Enhanced Resolution of Shocks (MOVERS) de¬ 
veloped by Jaisankar and Raghurama Rao Q is capable of capturing steady discontinuities like shock waves and 
contact discontinuities aligned with grid lines exactly without any numerical diffusion and is otherwise a low nu¬ 
merical diffusion algorithm. The exact discontinuity capturing, a feature shared by very few of the upwind schemes 
which made them popular, is a result of enforcing the Rankine-Hugoniot jump condition directly in the discretization 
process in MOVERS. With this feature the central discretization methods can compete in accuracy with the best of the 
Riemann solvers, in spite of being much simpler. However, the two versions of MOVERS presented by Jaisankar and 
Raghurama Rao ||T|, namely MOVERS-n (an n-wave model for the coefficent of numerical diffusion) and MOVERS-1 
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(a 1-wave model for the coefficient of numerical diffusion) require an entropy fix to avoid non-smoothness in the ex¬ 
pansion regions. To overcome the problem of an entropy fix, another version, MOVERS-L, is presented in that work 
based on the idea of shifting to base-line Rusanov solver Gl or Local Lax-Friedrichs (LLF) method 18] |9l in smooth 
regions and to MOVERS-n or MOVERS-1 near large gradients based on a limiter function. It is worth noting that 
while MOVERS-L overcomes the problem of using an entropy fix, due to the base-line solver being LLF method, a 
high dose of numerical diffusion will be present in the smooth regions while still capturing the discontinuities very 
accurately. In this work, we propose an alternative method to reduce the numerical diffusion in the smooth regions 
in MOVERS-L, based on appropriate guidance from the discretization of the entropy conservation equation. For a 
detailed mathematical analysis of entropy stability the reader is referred to the research work of Tadmor no). In this 
work, however, we present a simpler way of introducing entropy stability based on an LLF type discretization of the 
entropy conservation equation and using it for guidance while fixing the numerical diffusion in smooth regions for 
MOVERS-L, thereby introducing a new version of the algorithm which we name as MOVERS-LE. 

2. Numerical Method 

Before presenting the new entropy stable version of the algorithm, MOVERS-LE, its base-line solvers, MOVERS 
and MOVERS-L, are first introduced briefly in the following. 

2.7. MOVERS for 1-D Euler equations 

Consider the 1 -D Euler equations given by 


dU dF 
dt dx 


= 0 


where the conserved variable vector U and the inviscid flux vector F are defined as 


( 1 ) 
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Here, p is the density, u is the velocity, E is the total energy, given by £ = piy-\) T P pressure. The flux 
Jacobian matrix defined by A = has three real and distinct eigenvalues, Ai - u-a,A2 = u and Aj - u + a. Applying 
a semi-discrete cell-centered finite volume method to this system of conservation laws gives the update formula 
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where the superscript n indicates the time level and the subscript j represents the grid points at the cell centres, with 
j±^ representing the cell interfaces. The cell interface flux is written as an average flux plus a diffusive flux given by 


FU-UFj + Fj,,)-^-aj,.(Ujy,-Uj) 


(3) 


where represents the coefficient of numerical diffusion. For a stable scheme the general interface flux is given by 

F,^^(F2_ + FR)-^a,AU (4) 

where L and R represent the left and right states of the interface, the first term on the right hand side (average flux) 
represents the central discretization of the flux terms and the second term is the diffusive flux, with aj being the 
coefficient of numerical diffusion at the interface. The intention here is to obtain a using R-H condition to resolve 
steady discontinuities exactly. In this setting the relevant physical speed is the speed of the discontinuity. This speed, 
s, is given by Rankine-Hugoniot jump condition. Here R-H condition is enforced at the cell interfaces to get the 
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Figure 1: Moving Shock 


relevant expression for coefficient of numerical diffusion. For a discontinuous wave with a left state and right state 
given by Ul and t/«, the Rankine-Hugoniot jump condition is given by 


AF = sAU 

where, AF - Fr- Fi (5) 

AU^Ur- Ul 


In the finite volume method, the variation of the conserved variable is assumed as a piecewise polynomial function 
within each finite volume. Hence the discontinuity in the conserved variable can occur at a cell interface. Let us apply 
the above R-H condition at a cell interface, j+\. 

(6) 

To obtain the interface flux, we split the above R-H condition into two parts. Let us consider a discontinuity between 
j and j + \ moving with a speed as shown in Fig Q. This discontinuity can move to either left or right. The 
jump condition does not indicate the direction of the movement of the discontinuity. To distinguish the positive and 
negative speeds, we can split the speed of the discontinuity into two parts, positive and negative. By this, left and right 
moving discontinuities are separated. We can now split the R-H condition also into two parts as 


(7) 

and 

(8) 

Note that the right hand side is not changed as in the above equation of three variables, namely s, U and F, only two 
variables can be varied. Subtracting ([^ from (0, we obtain 
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( 10 ) 


Now we get the interface flux in-terms of the speed of the discontinuity which is obtained using R-H condition as 


j+ 




( 11 ) 


3 












( 12 ) 


Comparing both the interface fluxes ([^ and (11 1 , we obtain 


Thus the coefficient of numerical diffusion is equal to the absolute value of the speed of the discontinuity. We can 
evaluate the value of Sj^i numerically which will be described in the following. 

Note that R-H condition is used to obtain coefficient of numerical diffusion in MOVERS. Thus, this scheme 
captures steady discontinuities (shock waves, contact discontinuities, slip surfaces) exactly. To ensure that the speed 
of the discontinuity falls within the range of the eigen-spectrum, wave-speed correction is devised, which is discussed 
in the paper IT] . 

R-H condition for a system of hyperbolic conservation laws is given by 


AF = sAU 


(13) 


where 


AF ^Fr- Fl and AU ^ Ur- Ul 


(14) 


Here,for 1-D Euler equations, AF and AU are 3x1 matrices, hence we take Sj^i to be a 3 x 3 matrix. We choose 
to be a simple diagonal matrix as follows. 


Thus we obtain 
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Note that the coefficient of numerical diffusion is to be fixed as the absolute value of the speed of the discontinuity 
as we derived and we thus now have three speeds of discontinuities, each for one of the conservation laws. Since 
we have n number of wave speeds representing n conservation laws as the coefficients of numerical diffusion, this 
algorithm is termed as MOVERS-n. This can be easily extended to any hyperbolic system of conservation laws. In 
smooth regions of the flow, AU =» 0 and the evaluated diffusion values can be unphysical. To avoid this numerical 
overflow, the wave speed correction strategy is devised. The wave-speed correction is done such that the coefficient of 
numerical diffusion is limited within the range of the eigen-spectrum of the flux Jacobian matrix. In this framework 
the wave-speed correction is slightly modified. When the left flux is equal to right flux and there is jump in conserved 
variables, then this case refers to steady discontinuity, for which the speed of the discontinuity is zero and hence the 
steady shock waves and contact discontinuities are captured exactly without any numerical diffusion. Considering all 
the possibilities discussed, the implementation of coefficient of numerical diffusion is discussed below. 
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where e is a small parameter. Wave speed correction is enforced as follows. 


(17) 


(18) 
(19) 

Here, and A,„i„ are maximum and minimum eigenvalues computed across the interface and are given as 


^ jx- 1 — ^max if ^ j-i- i ^ ^max 
^ j +1 '^min 11 ^ j+ i ^ '^min 


'^max — max{max{\u -I- a\, \u\, \u — aDi, maxQu + a\, \u\, \u — «!)«) 
dmin = max{min(\u + a\, \u\, \u — a\)i, min(\u + a\, \u\, \u — «!)*) 
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The MOVERS is tested on steady contact and steady shock test cases to demonstrate its ability to capture them 
exactly. Details of these test cases are explained in the Section [TTTT As shown in the Fig ([^), MOVERS captures 
steady the contact discontinuity exactly. The steady shock wave problem is taken from mu, as shown in the Fig ([^) 
and MOVERS captures the steady shock also exactly without any numerical diffusion. Modified shock tube problem 
is one of the test problems used to assess the numerical method for its entropy stability property. The details of the 
problem are given in the section [TL2 MOVERS when tested on this problem, gives a sonic glitch in the expansion 
region as shown in the Fig (|^). MOVERS being an accurate scheme, has less numerical diffusion, and thus requires 
an entropy fix to avoid non-smoothness in the expansion regions. The details of the entropy fix are explained in the 
paper m. It is worth noting here that Roe scheme IfT^ gives an unphysical expansion shock, whereas MOVERS 
produces a milder sonic glitch. By using an entropy fix the problem of sonic glitch is eliminated, which is evident 
from Fig (|^). With the use of an entropy fix, although sonic glitch problem is overcome, the ability of MOVERS 
to capture steady discontinuities exactly is lost. As demonstrated in the Fig. ([^ the steady discontinuities are then 
diffused. Motivated by the desire to overcome this problem, Jaisankar and Raghurama Rao HI proposed a hybrid 
scheme named MOVERS-L, which uses the diffusion based on Rusanov scheme 121 in smooth regions of the flow and 
MOVERS at the discontinuities, based on a limiter function which is discussed in the next section. 


2.2. MOVERS-L 

The motivation for this version of the MOVERS is to retain the property of exact discontinuity capturing in 
steady state while avoiding the entropy fix. The Rusanov (LEE) scheme [|2l is a simple central solver which satisfies 
entropy condition and avoids expansion shocks. The coefficient of numerical diffusion in the LEE method is fixed 
as the maxium of the eigenvalues of the flux Jacobian matrix, from the left and rights states of the cell interface. In 
MOVERS-E, the numerical diffusion based on EEF scheme is used in the smooth regions of the flow and the numerical 
diffusion based on R-H condition as in MOVERS is used at the large gradients. The switching over between smooth 
regions and large gradients of the flow is done using a minmod limiter IIT3l . The coefficient of numerical diffusion in 
this method, at the cell interface, is given below. 


- a movers - <Pir)(a movers - oille) ( 22 ) 

where Qmovers is the diffusion given by the expressions to and aiif is the maximum spectral radius of the 
flux Jacobian matrix from the left and right states of the cell interface, (pir) is the minmod flux limiter which acts as a 
switch and distinguishes between smooth regions and high gradients of the flow. 


(r^+i) = minmod{\,r],rj^^) 


(23) 
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In the above expression of r|, as the denominator is the difference in conserved variables, when t//+i =» Uj the value 
of the expression can lead to large values of rj. To avoid the numerical overflow, the following definition of is 
used. 


rj = signiUj^i - if \Uj^i - Uj\ < 6 

where S is small number. The multivariable minmod function is defined as follows 11. 


(25) 


minmod{x\, X 2 , ■ ■.) 


minj if xy > 0 V j 
' maxj |xj| if Xy < 0 V j 


(26) 


0 otherwise 

As expected, MOVERS-E captures steady discontinuities exactly which is shown in the Fig. The results are 
also compared with EEF and MOVERS with entropy fix. This scheme is tested on modified shock tube problem with 


5 












a sonic point in the expansion fan to test the entropy stability property, as shown in Fig (j^. As expected, MOVERS-L 
avoids generating an expansion shock. Although MOVERS-L captures steady discontinuities exactly and yet avoids 
using an entropy fix, it is diffusive in smooth regions as is evident from Eig ^ since we are using LEE based diffusion 
in the smooth regions of the flow. Developing a numerical method which captures steady discontinuities exactly, is 
accurate in smooth regions and yet is entropy stable is the motivation for the present work. 

2.3. Entropy conservation equation based numerical diffusion 

The main focus of this work is to enforce entropy stability in smooth regions, retaining sharp resolution of discon¬ 
tinuities. The expansion regions with sonic points require finite numerical diffusion to avoid any expansion shocks, 
while shocks need less diffusion to resolve them accurately. Ensuring both these conflicting requirements in a numer¬ 
ical method is a non-trivial task. The basic idea here is to use the entropy conservation equation as a guide-line to fix 
the coefficient of numerical viscosity. Considering the second law of thermodynamics together with the conservation 
of mass, momentum and energy equations, we can rewrite the 1-D Euler equations using an entropy formulation as 
follows imiBiiisi. 


dp d (pu) 
dt dx 

d(pu) d[p+pu^) 

dt dx 

d (pS ) d (puS ) 
dt dx 

where entropy S is defined by 

S = Cy log -I-conifanf (30) 

Note that the energy conservation equation is replaced by an entropy conservation equation after a little bit of algebraic 
manipulation of the equations. This formulation is valid only for smooth regions, without discontinuities. Eor flows 
with discontinuities, the entropy conservation equation has to be replaced by an inequality as follows. 


(27) 

(28) 
(29) 


d{pS) dipuS) 


dt 


dx 


> 0 


(31) 


To keep the entropy bounded, mathematicians define the entropy with a negative sign so that we obtain the following 
equation. 

(32) 


dipS) ^ d(puS) ^ ^ 


dt dx 

In the absence of discontinuities (for smooth regions of the flow), the equality holds and the entropy conservation 
equation is given as 

d(n:^^ dtnut^^ 

(33) 


d(pS) ^ d(puS) ^ ^ 


dt dx 

Tadmor lfTOl has analyzed entropy stability and entropy conservation properties of numerical methods. He argues 


that, for any scheme to be entropy stable, it should satisfy additional inequality (32i. The basic strategy of the 
present research work is to use the above entropy conservation equation to get some information about how to fix 
the coefficient of numerical diffusion in finite volume method in smooth regions of the flow. Our approach is much 
simpler than that of Tadmor, while Tadmor’s approach is mathematically rigorous. 

A very efficient and yet simple strategy is to apply the LEE method for the entropy conservation equation to obtain 
the numerical diffusion and to simply use the same numerical diffusion for discretizing the Euler equations, using a 
finite volume method. Since the entropy conservation equation is applicable for smooth regions and moreover it is 
used in place of the energy conservation equation, the hope is that discretization will yield the required amount of 
numerical diffusion to yield a stable scheme for all the conservation equations in the smooth regions of the flow. It 
is also expected to yield an entropy stable scheme for Euler equations. Eventually, the plan is to couple this strategy 
to the strategy of MOVERS being used near the large gradients, thereby retaining the exact disctoninuity capturing 
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in a simple central discretization framework. Since the numerical diffusion in the LLF method is simply the local 
maximum of the eigenvalues of the flux Jacobian matrix for any system of equations, we first determine the wave- 
speed for the entropy conservation equation, which is given by 


Oes 


dipuS) 

d{pS) 


(34) 


Thus, the numerical diffusion rquired for applying the LLF method to the entropy conservation equation will be 


aiLF-ES = max[\uE\, |mr|] 


(35) 


In ges, the subscript ES refers to the phrase Entropy Stable (ES), where L and R denote the left and right states of 
a cell interface. Since this numerical diffusion is from a stable discretization of the entropy conservation equation 
for smooth regions, we can use the same numerical diffusion for discretizing Euler equations, with the expectation of 
obtaining a stable numerical method for smooth regions. This method is termed as MOVERS-E. This scheme is tested 
on steady shock problem and shock tube problem with sonic point. As shown in the Eig (|^) the scheme has no sonic 
point problem, and yet is accurate. But the scheme diffuses the steady shocks Eig (j^). By following the strategy of 
MOVERS-L, we device a hybrid scheme named MOVERS-LE (L refers to limiter and E refers to entropy stability) 
as described in the next section. 


2.4. Entropy stable hybrid solver 

The special discretization described in section 2.3 is an accurate and stable numerical scheme for only smooth 
flows as it diffuses the discontinuities. Near discontinuities, we can use the numerical diffusion based on Rankine- 
Hugoniot condition, as in MOVERS for sharp resolution of the high gradients. This hybrid scheme, named MOVERS- 
LE, can be obtained by utilizing the following coefficient of numerical diffusion at the interface. 


^MOVERS-LE,j+^ - ^MOVERS,i+\ + ^(O+P 


LLE-ES,j+{ 


■ ^ MOVERS. 


(36) 


where, aiiF-ES given by (35i and a movers is the R-H condition based numerical diffusion, as explained in 
Section 2.2 The results of this hyBrid scheme, MOVERS-LE, is presented in the next section for various test cases. 
In two dimensional case, a finite volume method is used with the above strategy being implemented across any cell 
interface in the direction normal to the interface. 


3. Results and discussions 

3.1. ID-Euler equation results 

MOVERS-LE is tested on following ID-Euler test cases, which are taken from Toro Oil to assess the numerical 
method for accuracy and robustness. 

3.1.1. Steady shock and contact discontinuity 

Eor the steady contact discontinuity, the initial condition are given as 

(p,M,p)fi = (L4,0.0,1.0) and (p,M,p)z, = (1.0,0.0,1.0) 

The initial discontinuity is at Xg - 0.5 and the final computation time is at f = 2.0. As shown the Eig(|^), MOVERS- 
LE captures this exactly. The steady shock problem is taken from Zhang and Shu IfTTI . This test case is built using 
R-H conditions. MOVERS-LE captures steady shock also exactly without any numerical diffusion, as shown in the 
Eig 

3.1.2. Sod’s modified Shock tube 

This test case is the modified version of Sod’s CD test case. The solution of this problem consists of shock and 
contact discontinuities which are moving to the right and an expansion with sonic point to the left. This test problem 
is used to assess the numerical method for its entropy stability. As shown in the Eig(|^), Roe scheme ifT^ gives an 
unphysical expansion shock at the sonic point, whereas MOVERS-LE gives a smooth expansion profile. 
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Figure 2: Steady Contact Discontituity (a) and Steady Shock (b) results with MOVERS using 100 grid points 
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Figure 3; Steady Contact Discontinuity (a) and Steady Shock (b) results with MOVERS and an entropy fix using 100 
grid points 


3.2. Strong shock 

This test case is used to check the robustness of the numerical method as it involves a strong shock wave of shock 
Mach number 198. The solution of this problem has a contact discontinuity, a strong shock and an expansion wave. 
The Fig. (|^) shows the solution of this problem using the proposed scheme, which demonstrates its robustness. 

3.3. Strong discontinuities 

This is severe test case as the solution of the problem has three discontinuities. This test case is also used to test 
the numerical method for accuracy and robustness. As shown in the Fig([T0^), MOVERS-LE does not encounter any 
problem and is robust. 
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Figure 4: First order result with MOVERS: (a) without and (b) with entropy fix for modified shock tube problem with 
sonic point using 100 grid points 
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Figure 5; Comparison of MOVERS-L results with MOVERS and LEE for Steady Contact Discontinuity (a) and 
Steady Shock (b) 


3.4. Slowly moving contact 

This test case is designed to test the numerical algorithm for its ability to resolve slowly moving contact dis¬ 
continuities. The solution of this test problem has a shock wave and a slowly moving contact discontinuity to right 
and a rarefaction wave to the left. As shown in the FigtfTO]^), the MOVERS-LE captures the slowly moving contact 
discontinuity quite accurately and its accuracy is comparable to that or Roe scheme. 

3.5. Test cases for 2D-Euler equations 

The present scheme, MOVERS-LE, is tested on various benchmark problems for 2D-Euler equations. The prob¬ 
lems are chosen to assess the numerical scheme for accuracy, robustness and ability to avoid shock instabilities. 
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Figure 6: Comparison of MOVERS-LE result with LEE and MOVERS (with entropy fix) for modified shock tube 
problem with sonic point using 100 grid points 



Eigure 7: Comparison of MOVERS-E results with MOVERS and LEE: (a) Steady Contact Discontinuity, (b) Modified 
shock tube with sonic point 
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Figure 8: Comparison of MOVERS-LE results with LEE for: (a) Steady Contact Discontinuity, (b) Steady Shock 
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Figure 9: (a) Comparison of MOVERS-LE result with LEE and Roe schemes for modified shock tube problem with 
sonic point using 100 grid points, (b) Strong Shock 


3.5.1. Grid aligned slip flow 

In this test case lfTOl . a Mach 3 flow slips over a Mach 2 flow, with the pressure gradient being zero and no jump in 
density. This slip flow problem tests the schemes for its ability to exactly capture grid aligned discontinuities. The left 
side of the domain is prescribed with inflow values and at all other boundaries, values are extrapolated from inside. 
Eig. (Ill shows the results with the present scheme, which exactly captures the slip stream as compared with LEE 
method. 


3.5.2. Oblique shock reflection 

In this test case Il20l an oblique shock is introduced at the top left corner of the computational domain. This 
oblique shock hits the solid wall and reflects back. The initial conditions for this problem are as given below. 

(P,m,v,p)o,,,, = (1.0,2.9,0,1/1.4) 
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Figure 10; (a) MOVERS-LE for Strong Discontinuities (b) Comparison of MOVERS-LE and Roe scheme for Slowly 
Moving Contact Discontinuity 



Eigure 11; Comparison of grid aligned slip flow problem results for MOVERS-LE (right) with LEE method (left) on 
a grid with 40 x 40 grid points: Mach contours (2.0; 0.05 ;3.0) 


(p,M,v,p)^_i_, = (1.69997,2.61934, -0.50633,1.52819) 

At the left and top sides of the domain inflow and post shock boundary conditions are prescribed as given above. 
Solid-wall and supersonic outflow boundary conditions are prescribed at the bottom and right boundaries. The results 
for the oblique shock reflection with MOVERS-LE for both first and second order accuracy are shown in the Eig. 
and Eig. ( [T3] ) respectively. As can be seen from the figures, the present scheme is slightly more accurate compared to 
MOVERS, and the contours are smoother. 
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3.5.3. Ramp in a channel 


In this test case ED a Mach 2 flow encounters a wind tunnel with a 15® ramp. The flow feature of this problem 
includes an oblique shock which reflects from the upper wall and interacts with the the expansion fan generated at 
the corner of the ramp. Fig. ( [T4| and Fig. (15 1 show the first and second order accurate results of MOVERS-LE 
compared with MOVERS. Although there is no noticeable difference in accuracy, it is worth noting that the contours 
are smoother with the use of present scheme. 



Eigure 14; Comparison of first order results of MOVERS (left) with MOVERS-LE (right): pressure contours 
(0.7:0.1:2.9) for the grids: (a) 60x20, (b) 120x40 and (c) 240x80 


3.5.4. Hypersonic flow past a half cylinder 

This is one of the multidimensional problems used to test the numerical algorithm for shock instability. Many 
upwind schemes generate the infamous carbuncle shock as demonstrated by Quirk El and Meng-Sing Liou Il22l . As 
shown in the Eig. ( [T6| ), the present method does not encounter any problem in capturing the flow features. 

5.5.5. Forward facing step 

In this unsteady test case Il2^ . a Mach 3 flow enters the wind tunnel with a forward-facing step. At time t - 4.0 
units, a lambda shock is visible with a triple point from which a slip surface is seen. The reflected oblique shock from 
top wall interacts with the strong expansion generated at the corner of the step. The left side of the computational 
domain is prescribed with inflow conditions. Supersonic outflow boundary condition is enforced at the exit and wall 
boundary conditions are prescribed at top, bottom and on the step. Eig. shows the first and second order accurate 
results for this problem with MOVERS-LE. As shown in the figure, the present scheme captures slip stream even with 
the first order accuracy, whereas many diffusive scheme like LEE do not capture it with first order accuracy. 

3.5.6. Shock diffraction 

This is another test case Il24l which assesses the scheme for shock instability. This test case has complex flow 
features involving planar moving shock moving with incident Mach number, a diffracted shock wave around the corner 
and a strong expansion wave. This strong expansion wave accelerates the flow and interacts with post shock fluid to 
further complicate the flow. Other distinct flow features are a slipstream and a contact surface. Numerical studies for 
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Figure 15: Comparison of second order results of MOVERS (left) with MOVERS-LE (right): pressure contours 
(0.7:0.1:2.9) for the grids: (a) 60x20, (b) 120x40 and (c) 240x80 



Eigure 16: Comparison of Eirst order results of MOVERS-LE with LEE on 45 x45 grids: density contours (2.0: 0.2: 
5.0) (a) LEE for Mach 6 flow (b) MOVERS-LE for Mach 6 flow (c) LEE for Mach 20 flow (d) MOVERS-LE Mach 
20 flow 


this problem are given by Hillier ll25l . Godunov-type and Roe schemes are known to fail for this case |l2j, as they 
produce oscillations at the planar shock. As can be seen from the Eig. (18l, the present method does not encounter 
any shock instability. It can also be noted that slipstream is seen with first order simulation itself, whereas second 
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Figure 17: Density contours (1.0: 6.5: 0.15) for forward facing step problem on 240 x 80 grids with MOVERS-LE 
for both first order (left) and second order (right) accuracy 

order results also resolve contact surfaces. 



Eigure 18: Density contours (0.5: 7.0: 0.25) for shock diffraction problem on 400 x 400 grids with MOVERS-LE for 
both first order (left) and second order (right) accuracy 


4. Conclusion 

A entropy-stable, hybrid numerical method for inviscid Euler equations is presented. This algorithm being a cen¬ 
tral solver, is free of complicated Riemann solvers and is not heavilty dependent on eigen-structure unlike upwind 
methods. This scheme uses the R-H condition-based numerical diffusion at the discontinuities and entropy conser¬ 
vation equation-based numerical diffusion in the smooth regions of the flow to capture grid aligned steady contact 
discontinuities and shock waves exactly, which is demonstrated through various benchmark problems. The idea of 
using numerical diffusion based on entropy conservation equation yields an entropy stable scheme for the smooth re¬ 
gions which is demonstrated using benchmark problems. This hybrid scheme is also tested on strong shock test cases 
to prove its robustness and its accuracy is comparable to the best of the approximate Riemann solvers. This entropy 
stable new algorithm is not prone to shock instabilities such as carbuncle shocks and avoids expansion shocks. This 
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new algorithm, MOVERS-LE, captures steady contact discontinuities and shock waves exactly, which is one of the 
desired features for a scheme in simulating viscous flows, an aspect which is currently being pursued. 
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